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Abstract 

Context. Most upcoming CMB polarization experiments will use direct imaging to search for the primordial gravitational waves 
through the B-modes. Bolometric interferometry is an appealing alternative to direct imaging that combines the advantages of inter- 
ferometry in terms of systematic effects handling and those of bolometric detectors in terms of sensitivity. 

Aims. We calculate the signal from a bolometric interferometer in order to investigate its sensitivity to the Stokes parameters paying 
particular attention to the choice of the phase -shifting scheme applied to the input channels in order to modulate the signal. 
Methods. The signal is expressed as a linear combination of the Stokes parameter visibilities whose coefficients are functions of the 
phase-shifts. 

Results. We show that the signal to noise ratio on the reconstructed visibilities can be maximized provided the fact that the phase- 
shifting scheme is chosen in a particular way called "coherent summation of equivalent baselines". As a result, a bolometric interfer- 
ometer is competitive with an imager having the same number of horns, but only if the coherent summation of equivalent baselines 
is performed. We confirm our calculations using a Monte-Carlo simulation. We also discuss the impact of the uncertainties on the 
relative calibration between bolometers and propose a way to avoid this systematic effect. 
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Introduction 

Measuring precisely the polarization of the Cosmic Microwave 
Background (CMB) is one of the major challenges of contempo- 
rary observationnal cosmology. It has already led to spectacular 
results concerning the cosmological model [Kovac et al., 2002 
Readhead et al. 2004| |Dunkley et al., 2008|lNoita et al., 2 008 
Ade et al., 2008) describing our Universe. Even more challeng- 
ing is the detection of the so-called B-modes in the CMB polar- 
ization, associated with pure tensor modes originating from pri- 
mordial gravitational waves enhanced by inflation. Discovering 
these modes would give direct information on inflation as 
the amplitude of the B-modes is proportional to the tensor 
to scalar ratio for the amplitude of the primordial density 
perturbations which is a direct product of inflationary sce- 
narii ]Liddle and Ly th, 2000 1 . Furthermore, it seems that most 
of the inflationary models arising in the context of string theory 
(brane inflation, ...) predict an undetectably small scalar to ten- 
sor ratio [Kall osh and Linde, 20071 . The discovery of B-modes 
in the CMB may therefore appear as the only present way to 
falsify string theories. Cosmic strings and other topological de- 
fects are also sources of density perturbations of both scalar and 
tensor nature. They are however largely dominated by the adi- 
abatic inflationary perturbations in TT, TE and EE power spec- 
tra and therefore hard to detect. It is only in the B-mode sector 
(BB power spectrum) that the tensor topological defects pertur- 
bation could be large [Bevis ~eTal., 2007| and have a different 
shape | Urrestil la et al., 2008| from those originating from infla- 
tion and hence be detectable |Pog osian and Wyman, 2007] |. 

Unfortunately, the inflationary tensor to scalar ratio seems 
to be rather small so that the B-modes are expected at a low 
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level as compared to the E-modes. The quest for the B-modes 
is a therefore tremendous experimental challenge: one requires 
exquisitely sensitive detectors with an unprecedented control 
of the instrumental systematics, observing at a number of dif- 
ferent frequencies to be able to remove foreground contami- 
nation. Various teams have decided to join the quest, most of 
them with instrumental designs based on the imager concept 
(BICEP, EBEX, QUIET, SPIDER, CLOVER). Another possible 
instrumental concept is an interferometer that has many advan- 
tages from the point of view of systematic effects (no optics for 
instance) and that directly measures the Fourier modes of the 
sky. Let us recall that the first detections of polarization of the 
CMB were performed with interferometers ]Kovac et al., 2 002 
Readhead et al., 2004] . Interferometers are however often con- 
sidered as less sensitive than imagers mainly because of the addi- 
tional noise induced by the amplifiers required for heterodyne in- 
terferometry whereas imagers use background limited bolome- 
ters. Another drawback of heterodyne interferometry is that it 
requires a number of correlators that scales as the square of the 
number of input channels limiting the number of channels actu- 
ally achievable | |CMB Task Force report, 2006) . 

A new concept of instrument called "Bolometric 
Interferometer" is currently under developpement 
(MBI HTimbieetal., 2003) , BRAIN (Polenta et al., 2003| 
Charlassie r et al., 2008| ). In such an instrument, the interference 
fringes are "imaged" using bolometers. We believe that such 
an instrument could combine the advantages of interferometry 
in terms of systematic effects and data analysis and those of 
bolometers in terms of sensitivity. The goal of this article is 
to investigate ways to reconstruct the Fourier modes on the 
sky (the so-called visibilities) of the Stokes parameters with a 
bolometric interferometer. In particular, we focus our attention 
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on the necessary phase-shifting schemes required to modulate 
the fringe patterns observed with the bolometer array. We show 
that one can construct phase-sequences that allow to achieve 
an excellent sensitivity on the visibilitites: scaling as y/Nh/N eq 
(where Nf, is the number of horns and N eq is the number of 
couples of horns separated by identical vectors hereafter called 
equivalent baselines) whereas it would scale as y/Nj,/ 
a non optimal phase-shifting sequence. 

This article is organised as follows: in section[T]we describe 
the assumptions that we make on the hardware design and on 
the properties of the various parts of the detector. In section [2] 
we describe how the signal measured by such an instrument can 
be expressed in terms of the Stokes parameter visibilities. We 
show how to invert the problem in an optimal way in section [3] 
and show how the phase-shifting scheme can be chosen so that 
the reconstruction is indeed optimal in section [4] We have val- 
idated the method we propose using a Monte-Carlo simulation 
described in section [5] We end up by some considerations about 
systematic effects induced by cross-calibration errors and pro- 
pose a way to avoid them in section|6] 



1. Bolometric Interferometer design 

In this section we will describe the basic design we assume for 
the bolometric interferometer and how the incoming radiation 
is transmitted through all of its elements. This will lead us to a 
model of the signal that is actually detected at the output of the 
interferometer. A schematic view of the bolometric interferome- 
ter is shown in Fig. [TJ 

1.1. Horns 

We assume that we are dealing with an instrument which is ob- 
serving the sky through Nh input horns placed on an array at 
positions d,. All horns are supposed to be coplanar and looking 
towards the same direction on the sky. They are characterized by 
their beam pattern on the sky noted B\ n (n) where n is the unit 
vector on the sphere. Two horns i and j form a baseline which 
we label by < b < Nh(Nh - 1)12- 1. The phase difference be- 
tween the signal reaching the two horns from the same direction 
n of the sky is such that: 

Ej(n) = Ej(n) exp(2 iniib ■ n), where Ub = (dj - df)/A, (1) 
where A is the central observing wavelength. 

1.2. Equivalent baselines 

It is clear that if two baselines b and b' are such that Ub = uy, 
then the phase shifts associated with the two baselines are equal, 
a fact that we shall extensively use in the following. All base- 
lines b such that Ub = Up form a class of equivalent baselines 
associated with mode Up in visibility space. For all baselines b 
belonging to the same class B, the phase difference between the 
two horns i and j is the same: 



Ej(n) = Ei(n)e\p(2inup ■ n). 



(2) 



The number N± of different classes of equivalent baselines de- 
pends on the array, and the number of different baselines in an 
equivalence class also depends on the particular class. For in- 
stance, if we consider a square array with Nh = A^ ide horns, there 
are N± = 2 A^ s id e (A^ s id e - 1) classes, and the number of equivalent 
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Figure 1. Schematic view of the bolometric interferometer design con- 
sidered in this article. 



baselines in the class associated witlQ 

1 I with 1 < I < N, ide - 1 for m = 
m j 

and -|N side - 1| < / < N side - 1 for 1 < m < N side - 1, 
is Neq(J3) = (N slde - IflXAfcde - m). 



1.3. Polarization splitters 

In order to be sensitive to the polarization of the incoming 
radiation, we also assume that at the output of each hom 
there is a device which separates the radiation into two or- 
thogonal components noted || and ±. Such a separation can 
be achieved with an OrthoMode Transducer (OMT) in wave- 
guide IPisano et al., 2007|, finline ||Chattopaddhyay et al. 1999| 
or planar |Engargiola et al. 1999 IGrimes et al. 200711 technolo- 



gies. Each horn therefore has two outputs measuring the electric 
field integrated through the beam in the two orthogonal direc- 
tions. The contribution coming from direction n for polarization 
77 (|| or ±) is: 



e? («) = B in (n)EKn). 



1.4. Phase-shifters 



(3) 



Important components of the required setup are the phase- 
shifters placed on each of the outputs that allow the phase of the 
electric field to be shifted by a given angle that can be chosen and 
controlled externally. This is the way the signal is modulated in 
order to recover the incoming information. For now we do not 
make any assumptions on the possible values of the angles but 
we will see that they have to be chosen carefully in order to op- 
timize the signal to noise ratio. The signal after phase-shifting 
coming from direction n with polarization 77 is: 



e/'in) = cxp(i<f>1)e'/(n). 



(4) 



For obvious hardware reasons, all phase-shifters in the setup 
have to be identical and deliver the same possible phase-shifts. 

1.5. Beam combiner 

In order to be able to perform interferometry, the beam of each 
horn has to be combined with all the others so that all pos- 
sible baselines are formed. The realization of a beam com- 
biner is an issue in itself that will not be assessed in the 
present article. As an example, this can be achieved using 



in units of the smallest baseline in the array. 
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a Butler combiner [Da lTOmo, 20 03 1 or with a quasi-optical 
Fizeau combiner such as the one used for the MBI instru- 
ment |Timbie et al., 2003 1. All of these devices are such that the 
2Nh input channels result after passing through the beam com- 
biner in A^out output channels that are linear combinations of the 
input ones. To be able to conserve the input power in an ideal 
lossless device, the number of output channels iV out has to be at 
least equal to the number of input channels 2Nh- In the output 
channel k the electric field is: 



Zk(n) 




(5) 



1=0 J7=0 



where the yl, coefficient model the beam combiner, r\ — 1 or 
respectively corresponds to || and ± polarizations. We choose 
to deal with configurations where the incoming power is equally 
distributed among all output channels, meaning that the coef- 
ficients y l ki have a unit modulus: p^jCfc) = 1 . In order to sim- 
plify the notation, we include the y n k . phases in the phase-shifting 
terms as . = (p'j + Arg(y^ .) so that: 



Zk(n) = 



N h -l 1 
out (=0 ;;=() 



(6) 



1.6. Total power detector 

The signal from each of the outputs of the combiner is not de- 
tected in a coherent way as in a heterodyne interferometer but 
with a bolometer through its total power averaged on time scales 
given by the time constant of the detector (larger than the fre- 
quency of the EM wave). We assume that the bolometers are 
background limited, meaning that their noise variance is propor- 
tional to their input power. The power on a given bolometer is 
then: 

' ) (7) 

' time 



Zk(n)dn 



= / {z k (n)zt(n)\ ime dndn'. 



(8) 



The signal coming from different directions in the sky are inco- 
herent so that their time averaged correlation vanishes: 



(z k (n)zt(n'j) = (\z k (n)\ 2 ) 6(n - «') 

\ * / time \ / time 

= \z k (n)\ 2 5(n-n'). 



(9) 
(10) 



From now on, z is then implicitely replaced by its time-averaged 
value. The signal on the bolometers is finally: 



Sk = J \z k (n)\ 2 dn. 

2. Stokes parameter visibilities 



(ID 



Developping the signal on the bolometers in terms of the incom- 
ing electric fields easily shows autocorrelation terms for each 
channel as well as cross-correlation terms between all the possi- 
ble pairs of channels: 

2 



+2Re 



i<j mm 



Z<(«)ex P (^,) 

ij=0 

*(n)ef *(n)ex P (j(0^-0^.)) 



>d«.(12) 



The electric fields from different horns are related through Eq. [2] 
and introduce the Stokes parameters that are generally used to 
describe a polarized radiation: 

(|£J 2 ), (13) 

Q= (|£||| 2 )-(|£J 2 ), (14) 

U= (EnE*) + (E*E x ) =2Re(£||Z?t), (15) 

V = i((E } \E* ± ) - (e;e ± )) = -21m . (16) 

The Stokes parameter visibilities are defined as {S stands for /, 
Q, U or V): 



(17) 



Vs (up) = J S (n)B 2 n (n) exp{2imip • n)dn. 

The phase-shift differences for a baseline b formed by horns ; 
and j measured in the channel k are: 



AO- 


= *» . - O ., 

k,i k,j' 


(18) 




= Ki-Kp 


(19) 




= <i-Kp 


(20) 






(21) 


Putting 


l all these definitions into Eq. 12 and after 


some cal- 



culations one finds that the signal on the bolometer k can be 
expressed purely in terms of the Stokes parameter visibilities 
and the phase-shifting values (the subscript b stands for all the 
Nh(Nh - l)/2 available baselines and n k is the noise): 

tf»(%-l)/2-l 

S k = A k • S + r ^ ■ ^ + W 

b=0 

where the first term is the autocorrelations of all horns and the 
second one contains the cross-correlations, hence the interfer- 
ence patterns. We have used the following definitions: 



AL 



■z 



S' 



1 



cos AOj I + cos AO^ 



-(sin A<J>;l 
cos AO' 1 1 ■ 

k.b 

— (sin A<dJ|J - sinAflv 



+ sinAO^) 
- cos AO^ 

k,b ' 



cos AO"^ + cos A<b kb 



/ I(n)B 2 (n)dn 
J Q(n)B 2 (n)dn 
J U(n)B 2 (n)dn 
J V(n)B 2 (n)dn ) 



Re MOB*)] 
lm[y,(u b )] 
Re [V e (« 4 )] 
Im[V e (« fc )] 
Re lV v (u b )] 
Im [Vu(u b )] 
Re[Vv(«(,)] 
Im[V v («„)] 



(23) 



'Vt 



(24) 



-(sinAO^ + sinA*^) 
-(sinAO^-sinAO^) 
-(cos AcD^-cos A<S>fl 

All of this can be regrouped as a simple linear expression in- 
volving a vector with all the sky informations (Stokes parameter 
autocorrelations S and all visibilities 'V/,) labelled X and another 
involving the phase-shifting informations (A* and T k b) labelled 
A k : 



A k -X + n k . 



(25) 



Finally, various measurements of the signal coming from the dif- 
ferent channels and/or from different time samples with differ- 
ent phase-shifting configurations can be regrouped together by 
adding columns to A which then becomes a matrix A and trans- 
forming the individual measurement S k into a vector S: 

S = A-X + n. (26) 
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3. Reconstruction of the visibilities 

Once one has recorded enough data samples to invert the above 
linear problem (we will call such a period a sequence in the fol- 
lowing), the solution is the usual one assuming that the measure- 
ments noise covariance matrix is N - (n ■ n'): 



X = (A t -N~ i - Ay 1 - A' ■ N~ l S, 
with covariance matrix: 

N = l(x- (x)) ■ (x - (!))') = (A' ■ AT 1 . Ay 



(27) 



(28) 



3. 1 . Regrouping equivalent baselines 

One sees that the dimension N„ of the X vector of unknowns is 
rather large: N u = 3+SxN b where N b - N h (N h -l)/2 is the num- 
ber of baselines formed by the input horn array. For a large horn 
array this number can become really large. A 10 x 10 array has 
for instance N b = 4950 baselines and N„ = 39603 unknowns. 
One needs at least as many data samples as unknowns (and in 
many cases more than that) so this would involve manipulations 
of very large matrices. In fact as we said before, depending on 
the relative positions of the input horns, there may be a lot of 
equivalent baselines: different couples of horns separated by the 
same vector Up hence measuring exactly the same visibilities. 
It is clearly advantageous to regroup these equivalent baselines 
together in order to reduce the dimension of the system. As we 
will see below there is a huge extra-advantage to do it this way 
in terms of signal-to-noise ratio if one chooses the phase-shifters 
angles wisely. 

In the case where the input horn array is a square grid with 
size A^side = *4Ni, the number of different classes of equivalent 
baselines is N t = 2N side (N side - 1) = 2{N h - = 180 for a 

10x10 horn array, hence reducing the number of unknowns to 
1443 which is a huge improvement. It is obvious that all equiv- 
alent baselines measure the same visibilities and can therefore 
be regrouped together in the linear problem leading to the same 
solution as considering the equivalent baselines separately. One 



just has to reorder the terms in Eq. 22 as first a sum over all dif- 



ferent baselines /3 and then a sum over each of the baselines bp 
equivalent to ft coming on the output line k: 



ArS + 



/J=0 



'Vb + n k , 



(29) 



changing the T vector to: 



2 



b B =0 



cos AO"" + cos AO 

k.be 



k.b„ 



-(sin AO" -t-sinAO^) 
cos AO 1 , 1 1 - cos A<bj~j- 



*k,b, 



Lb, 



(30) 



-(sin AO -sinAO^) 
cosAO'lt -HcosAO^ 1 ' 

k,bp k,ba 

-(sinAO"'; + sinAOf|) 
-(sinAO^ -sinAO^I) 

v k,bfj k,b„> 

-(cosAO'h -cos AO;* 

y k.be k.be 



Let's recall that each column of the matrix A corresponds to 
phase-shifters configurations encoded in Ykj} (for all different 
baselines /3) and A^. 



3.2. Coherent summation of equivalent baselines 

We shall now investigate the noise covariance matrix for the re- 
constructed visibilities and how one can possibly optimize it. We 
will assume for simplicity that the noise is stationary and un- 
correlated from one data sample to another where we call data 
sample the output of one of the N out output channel during one 
of N, time samples. Therefore there are Nd = N oul x N, data 
samples. If the photon noise corresponding to one horn on one 
detector measured during one time sample is cr , then the noise 
covariance matrix of the measured data samples is: 



N = 



x n, 



(31) 



where 11 is the Nd x Nd identity matrix. Assuming that the time 
variation of the A matrix can be neglectecQ It can be trivially 
extended to a Nd x A^, matrix Af,,. p = A^p. 
In terms of this extended A matrix, the visibilities covariance 
matrix (see Eq.|28]i writes: 



N = 



N nut 



(a'-a)~ 



(32) 



We have regrouped all equivalent baselines together in A, 
each of its elements is therefore the sum on Af eq sines and cosines 



of the phase-shifting angles (as expressed in Eq. 30 1. We will 
assume here that the angles are chosen randomly and uniformly 
from a set of possible values between and 2n. Now there are 
two possibilities depending on the choice for the phase-shifting 
angles for all baselines equivalent to a given one: they can all 
be different or they can all be equal. We refer to this choice as 
incoherent or coherent summation of equivalent baselines: 

- Incoherent summation of equivalent baselines: each of the 
sum of the two sine/cosine functions of the uniformly dis- 
tributed angles has zero average and a variance 1. Each ele- 
ment of Ykfi is the sum of Af eq of these and the Central Limit 
Theorem states that it will have zero average and a variance 

— Coherent summation of equivalent baselines: then each ele- 
ments of is j7—N eq (/J) times the same angle contribution 
with variance 1 . The matrix elements ends up having a vari- 



ance 



Coming back to A 1 ■ A, the multiplication by the transpose will 
add together all the Nd different data samples. The off-diagonal 
elements will cancel out to zero because the angles are uncor- 
rected from one channel to another. The diagonal elements will 
however average to the variance of the elements in A multiplied 
by Nd- So finally, depending on the choice between incoherent 
or coherent summation of equivalent baselines, the visibility co- 
variance matrix will scale in a different manner: 



N 



1 



N oat N d iV„,08) 



°V*i> i N l 



for incoherent summation, 



(33) 



for coherent summation. 



The latter scaling is clearly more advantageous and optimises 
the reconstruction of the visibilities. In fact this result is quite 
obvious: if the phase-shifting angles for equivalent baselines are 



2 Time variation of the A matrix is of course a source of systematic 
errors and must be studied as such. 
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all different, the coefficients of the linear problem that one wants 
to invert will always be smaller than if the summation of equiva- 
lent baselines is performed coherently. The signal to noise ratio 
on the visibilities will therefore be optimal if one maximises the 
coefficients, which is obtained by choosing the coherent summa- 
tion. 



3.3. Comparison with classical interferometers and imagers 

The variance on the visibilities obtained above in the case of a 
coherent summation of equivalent baselines can be rewritten: 



H 20 H 30 



cr; 



N h 



N eq (ft) N,N eq (ftY 



that can be comparecQ to 
in |Hobson and Magueijo, 1996) which 



(34) 

formula (28) 
is the equivalent 
for heterodyne interferometry crilN,N eq (ft). We see that the only 
difference introduced by bolometric interferometry is the factor 
Ni,/N eq (ft). In average, the number of equivalent baselines is 

^V eq ) = (N h (N h - l)/2)/N* - N h /4, but is much larger for small 
baselines. The design of the instrument has to be such that the 
"interesting" baselines are very redundant leading to a Ni,/N eq 
closer to one. The resulting expression of the variance on 
visibilities for bolometric interferometry therefore only differs 
by this slightly larger than one factor with respect to heterodyne 
interferometry. The important point is that the value of cr for 
bolometric interferometry is typical of a bolometer (photon 
noise dominated) hence smaller than what can be achieved with 
HEMT amplifiers in a heterodyne interferometer. 

This result can be summarized as follows: a bolometric in- 
terferometer using coherent summation of equivalent baselines 
can achieve the sensitivity that would be obtained with an het- 
erodyne interferometer with the noise of a bolometric instrument 
(and without the complexity issues related to the large number 
of channels). Such an instrument would therefore be competitive 
with an imager that would have the same number of bolometers 
as we have input channels in our bolometric interferometer. A 
detailed study of the comparison between a bolometric interfer- 
ometer and an imager from the sensitivity point of view is in 
preparation [Charlassier et al., in preparation | . On the opposite, 
if the equivalent baselines are summed incoherently, it is obvi- 
ous that the sensitivity would be very poor due to the absence of 
the 1 /N eq additional factor. 

The next section shows how it is possible to choose the 
phase-shifting sequences in such a way that the prescription of 
coherent summation of equivalent baselines is enforced. 



4. Choice of the optimal phase-sequences 

One wants the phase-shifting scheme to be such that equivalent 
baselines have exactly the same sequence but that different base- 
lines have different phase-shifts so that they can be disentangled 
by the linear inversion corresponding to Eq. 29 Now let's see 



how to comply with this constraint of having equivalent base- 
lines correspond to identical phase differences. An important re- 
mark is that, as can be seen in Eq. [6] the phase-shift have two 
different origins: the phase-shifters themselves on the one hand 



3 The notations are different : t vm in [Hobson and Magueijo, 1996 1 
has to be replaced by our N t , their n vis is the number of equivalent base- 
lines yV cq . In our article crn corresponds to s£l s in their article as a noise 
equivalent power NEP has to be replaced by NET x Q when talking 
about noises in temperature unit rather than in power unit. 
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Figure 2. Choosing all the phase sequences from that of one horn and 
two phase sequence differences (represented in green). 

whose angles can be chosen to follow a given sequence and are 
the same for all output channels and on the other hand, the phase- 
shifts coming from the beam combiner. Each input will be la- 
belled by the horn number < i < — 1 and the polarization 
direction rj. Each output is labelled by its number 1 < k < N out ■ 
The phase-shift differences are therefore 



k®k,m,jrf - Win ~ 4>jrf) + (ftkJn - ftkJn') ■ 
phase-shifters beam combiner 



(35) 



a) Phase-shifter phase differences (unpolarized case): 
We assume that the horns are placed on a square array with size 
^side = V^Vfi as m Fig- [2] In this case, the position of all horns 
can be parametrized, in units of the minimum horn separation, 

as a vector rf, = [ 1 | where /,■ and m,- are integers running from 



to iV s ide - 1 such that i - 1/ + A^ide/w,-. In this case, we have seen 
that there are N± = 2(Nh - y/NJ,) different classes of equivalent 
baselines labelled Up. Forgetting about polarization, the phase 
sequences can be constructed from a vector of two independent 
random phase sequences h(t) and v(t) which separate the hori- 
zontal and vertical directions in the horn array: 



<i(t) = di ■ s(f) where s(t) 



h(t) 
v(t) 



(36) 



The phase shift difference associated with the baseline between 
horns ; and j is 



(di - dj) ■ s(t), 



(37) 



and it is clear that the phase shift difference sequences will be 
the same for all baselines such that di - dj = Up, where ft is one 
of the classes of equivalent baselines. Because the two random 
sequences h(f) and v(t) have been chosen independent, the phase 
sequences associated with two different baselines classes ft + ft' 
will be different. 



b) Separating polarizations: 

Looking at formula ((24]), it is clear that one will not be able to 
separate V/ and V G visibilities, unless one uses two independent 
vectors of sequences S\\(f) + Sj_(0- However, in this case Vjy and 
Vy are not measured with maximum accuracy because the phase 
shift differences 



"j J. 



= *[[(*) ■ d t - s x (f) ■ d> 



are not equal for two different but equivalent baselines, so that 
they do not add coherently. One is therefore led to use alternately 
two measuring modes: 
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1 . One mode where S\\(f) + s ± (t), where phase shifts differences 
read: 

4>i n ~ 4>jn = ■ (dj - dj). 

In this mode, V/ and Vg are measured with maximum accu- 
racy (noise reduction oc N^ q ), but Vy and Vy are only mea- 
sured with noise a reduction oc Af eq . 

2. One mode with ,$||(f) = s ± (t) = s(t). Then however, one can- 
not measure Vy because 

<f>i\\ - (f>j± = s(f) ■ (dj - dj) = <p i± - 4>j\\, 

therefore one must introduce two more sequences cy(f) + 
Cj_(f) (one of them may be zero) independent from one an- 
other and from s(t), such that (f> ir) - s(t) ■ di + c n (t). Then: 

<Pi\\-</>j± = s(f)-(Ui-Uj) + c ll (t)-c ± (f) 



whereas 



S(t)-(Ui-Uj) + C ± (t)-C\\(f) 



but 



<Pi\\ ~ </>j\\ = 4>i± ~ <Pj± = s (t) ■ ("i - «j)> 

which means that V/, and Vy are measured with maxi- 
mum accuracy (noise reduction oc N^ q ), but Vg is not mea- 
sured at all. 

We mentioned before that we need all phase-shifters to be 
identical: they all have to be able to produce the same phase- 
shifts (let us call this ensemble W) but in different order. If h n (t), 
v v (f) and c n {f) are sequences of elements belonging to *¥, then 
the phase for any horn also has to belong to *P, meaning that 
0oj(O = // hjAt) + nii v n (f) + c n {t) has to belong to W. As shown in 
Appendix [A] this requires to choose the n$ values of the phase- 
shifts regularly spaced between and 2n as: 



In 



(n = 0,. ..,ns - 1). 



The elementary sequences h v (i), v n ( f) and c n (t) are uniform ran- 
dom sample of N s values taken among the elements of 
They must be chosen independent from one another to make 
sure that unequivalent baselines do not share the same sequence 
of phase differences. 



c) Beam combiner phase difference: 

As was said before there are two main designs for the beam 
combiner: Butler combiner |DalFOmo, 2003 1 or quasi optical 
combiner [Tim bie et al., 2003| . Without going into details, let us 
say that identical phase shifts for equivalent baselines are natu- 
rally obtained for the quasi optical combiner, and are achieved 
through an adequate wiring for the Butler combiner. 



d) Summary and expected accuracy: 

Finally, in order to recover the visibilities keeping to the "co- 
herent summation of equivalent baselines" criterion, one only 
has to build phase sequences that successively follow modes 1 
and 2 on an equal footing, build the corresponding A matrix and 
solve the system. There is a price to pay: during the first half se- 
quence, V/ and Vg are measured with optimal accuracly but Vy 
and Vy are not, during the second half sequence, V/, Vy, and Vy 
are measured with optimal accuracy but Vg is not measured at 
all. We therefore expect the sensitivity on Vg, Vy and Vy to be 
down by roughly a factor of V2 with respect to the sensitivity on 
V/, although the sensitivity on Vy and Vy will be slightly better 
constrained than on Vg. 



5. Monte-Carlo simulations 

We have investigated what was discussed above using Monte- 
Carlo simulations. There are three approaches that have to be 
compared for the reconstruction of the Stokes parameter visibil- 
ities. 

- Considering all baselines independantely without regrouping 
the equivalent ones. We expect this method to have error bars 



scaling as 1 / 



The system to solve is large in that case. 



- Regrouping the equivalent baselines together but without any 
choice for the phase-shifts so that they don't add in a coher- 
ent way. We expect this method to be exactly equivalent to 
the previous one but with a reduced size of the matrices. 

- Following the strategy to regroup equivalent baselines and 
choose the phases so that they are coherently added. We ex- 
pect the error bars to scale as 1 /N eq and therefore be the most 
efficient. 

In each case, we have simulated random visibilities with Vg, Vy 
and Vy a hundred times lower than V/ as expected from the CMB 
and calculated the signal expected on the bolometers using the 
phase-shift values for the three above strategies. We then added 
Gaussian noise with a variance o-^ c = cr^Nh/Nt to the bolometer 
signal. In each case we have performed a large number of noise 
and phase-shift sequence realisations. For each realisation, we 
have stored the reconstructed and input visibilities and analysed 
the residuals distributions. We have investigated the three above 
strategies and also the behaviour of the third one (coherent sum- 
mation of equivalent baselines) with respect to the two free pa- 
rameters: the length of the phase-shift sequence before inverting 
the linear problem and the number of different phase-shift angles 
(regularly spaced between and 2n as shown in Appendix [A). 



(38) 5.1. Scaling with the number of equivalent baselines 



We show in Fig.|3]the scaling of the RMS residuals on the vis- 
ibilities as a function to the number of equivalent baselines. We 
have divided the RMS by cr -\fNfjN~, in order to isolate the ef- 
fects that are specific to bolometric interferometry and depend 
on the way equivalent baselines are summed (see Eq.[33|. We see 
that as expected the scaling is oc 1 /N eq if one solves the problem 
by maximizing the signal to noise ratio using our coherent sum- 



mation of equivalent baselines. The poor 1 / 



scaling is also 



observed when all baselines are considered separately or when 
the phase-shift angles are not choosen optimally. 

5.2. Scaling with the number of samples and number of 
different phases 

Let's now concentrate on the optimized strategy described 
above: coherent summation of equivalent baselines. We show in 
Fig. [4] the scaling of the RMS residuals on the visibilities with 
respect to the length of the sequence and the number of different 
phases achieved by the phase-shifters (as shown in Appendix [A] 
these have to be regularly spaced between and 2n). The RMS 

values have been divided by ff- ^jf. 

One observes (Fig. |4]left) that the linear problem is singular 
when the number of different phases is not sufficient. Varying 
the number of horns in the array led us to derive the general 
scaling 2 y/Ni, for the minimum number of phases. Increasing 
the number of possible angles does not improve the residuals. 
Concerning the length of the sequence (Fig. |4] right), one ob- 
serves that when it is slightly larger than the number of unknows 
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Figure 3. Relative RMS on visibility residuals for V; (left) and Vg (right). The behaviour for and VV would be exactly the same as for Vy. 
The RMS on the residuals has been divided by the expected cr x scaling for each strategy to exhibit only the dependence with the number 
of equivalent baselines. The data points were fitted with linear slopes in Log-Log to measure the power of the scaling. One sees that the strategy 
where equivalent baselines are summed in a coherent way leads to a much better scaling oc -L than the other strategies that both scale as oc * — . 



(N u = 3 + 8 x N± where is the number of different baselines, 
N± = 2(Nh - yfNh) for a square array) then the reconstruction 
of the visibilities is not optimal due to the lack of constraints. 
Optimality is progressively reached when integrating a larger 
number of samples before inverting the problem. A reasonnably 
optimal result is obtained when Nd — 4xN u . The expected y2 
difference between the accuracy on V/ and that on Vq, Vy and 
Vy (due to the fact that we have have to perform two successive 
phase-shifting schemes in order to measure all three polarized 
visibilities) is also confirmed by the simulation. 

6. How to proceed with a realistic instrument ? 

When dealing with a realistic instrument one has to account for 
systematic errors and uncertainty to choose the precise data anal- 
ysis strategy. We do not want to address the wide topic of sys- 
tematic effects with bolometric interferometry in this article (we 
refer the interested reader to [Bu nn, 2007| where systematic is- 
sues for interferometry are treated in a general way) but just want 
to stress one point that is specific to the method we propose here, 
related to intercalibration of the bolometers in the detector array. 
Inverting the linear problem in Eq. 29 is nothing but ex- 



pressing the Stokes parameter visibilities as linear combinations 
of the Nd signal measurements performed with different phase- 
shifting configurations. These measurements can be those of the 
N oul bolometers each in N, time samples. This is where inter- 
calibration issues have to be considered. Linear combinations of 
signals measured by different bolometers are extremely sensitive 
to errors in intercalibration and will induce leakage of intensity 
into the polarized Stokes parameters if it is not controlled up 
to an exquisite accuracy. So we claim that combining different 
bolometers in the reconstruction of the visibilities in a bolomet- 
ric interferometer such as the one we describe here is not a wise 
choice unless the bolometers array is very well intercalibrated 
(through precise flat-fielding). The solution we propose is to treat 
all the bolometers independantly, inverting the linear problem 
separately for each of them. This requires a lot of time samples 
for the phase-shift sequences but is safer from the point of view 
of systematics. As a realistic example, for a 10 x 10 elements 



square input array, the number of different baselines is 180 and 
the number of unknowns is 1443. An optimal reconstruction of 
the visibilities can therefore be achieved with ~6000 time sam- 
ples. The duration of the time samples is driven by both the time 
constant of the bolometers (very short with TES) and the speed 
achieved by the phase-shifter to switch from one phase to the 
other. A reasonnable duration for the time samples is about 10 
msec which would correspond to sequences lasting about one 
minute. It is likely that the cryogenic system of such a bolomet- 
ric interferometer would ensure a stable bath on the minute time 
scale so that the knee frequency of the bolometric signal would 
be smaller than 1 min -1 . In such a case, the noise can be consid- 
ered as white (diagonal covariance matrix) during each sequence 
and the inversion gets easily tractable even with 6000 samples 
vectors. We are currently performing fully realistic simulations 
including systematic effects, the results will be presented in a 
future publication. 

Conclusions 

We have investigated the way to reconstruct the Stokes param- 
eter visibilities from a bolometric interferometer. It turns out 
that all three complex Stokes parameter visibilities can be recon- 
structed with an accuracy that scales as the inverse of the number 
of equivalent baselines if one follows a simple prescription: all 
equivalent baselines have to be factorized together in a coherent 
way, meaning that the phase-shift differences have to be equal 
for equivalent baselines. We have proposed a simple way to con- 
struct such phase-shift sequences and tested it on a Monte-Carlo 
simulation. The simulation confirms that the scaling of the errors 
on the visibilitites is oc ~\fNj,/N eq if one follows our prescription 
but JNh /N eq otherwise. 

The main conclusion of this article is therefore that a bolo- 
metric interferometer is competitive with an imager having the 
same number of horns (instrumental noise on the power spec- 
trum oc 1 /Nh) but only with an appropriate choice of the phase- 
shift sequences (coherent summation of equivalent baselines). 

We also discussed the data analysis strategy and proposed 
a solution to the possible cross-calibration issues between the 
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Figure 4. Scaling of the RMS residuals (divided by 



*■) on the Stokes parameter visibilities with respect to the number of different phases 



achieved by the phase-shifters and the number of different phase configurations used for the analysis (length of the sequence). One sees on the 
left that for a number of horns of 64 one needs at least 12 or 13 different angles to be able to solve the linear problem. It is clear from the plot 
on the right that the longer the sequence, the better the residuals, but a plateau is rapidly reached when the number of samples is around 4 times 
the number of unknowns in the linear problem. One can also see the factor = V2 between the accuracy on the intensity and the polarized Stokes 
parameters due to the two-steps phase-shifting scheme that we have to perform to be able to reconstruct them all. 



different bolometers. Even though one has simultaneously 7V out 
measurements of the signal with different phase configuration, 
it might be preferable not to combine these measurements but 
to reconstruct the visibilities on each bolometer separately and 
combine the visibilities afterwards. Such a strategy would in- 
crease the length of the phase-shifting sequences, but in a rea- 
sonable (and tractable) way thanks to the intrinsic shortness of 
our proposed phase-shifting scheme. 
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Appendix A: Proof of the necessity of having 
regularly spaced phase-shift values 



When we use the phase-shift configurations of Eq. 36 the an- 
tenna with coordinates (i, j) will be phase-shifted by: 



cpijit) = ih(t) + jv(t) + c(t). 



(A.l) 



In practice we are only able to construct a limited number of 
different phase-shifters, and the phase-shift sequences h{f), v(t) 
and c(?) will be independent random sequences of phase-shifts 
taken from the same set 3> of n phase-shifts <p p . For all phase- 
shifts in Eq. A. 1 to belong to O, it is necessary that lxcf) p (modulo 



In) also belongs to Let us write the smallest non-zero element 
of O as: 



2n 
n + e 



n e IN, 



< e < 1. 



(A.2) 



(n + l)<7> m in (modulo 2n) should also belong to <J, but 



27Z"( 1 € ) 27Z"( 1 f ) 

(n + l)</>mi n -2n -\ — = — (modulo 2n). (A. 3) 



n + e 



n + £ 



Therefore (n + 1 )(f> n 



(modulo 2n), and cannot belong to 



O unless e = 0. One concludes that the set 0„ of n phase-shifts 
has to be of the form: 



<D„ 



>n,p 



2np 
n 



n e IN, peN, < p < n 



(A.4) 



which finally is a quite obvious choice. 



